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Abstract 

The process e + e~ — * fi + u~~ has been studied with SND detector at VEPP-2M e + e~ col- 
lider in the vicinity of </>(1020) resonance. The product of branching ratios of <f> meson into 
leptons 



y/B{<j> ->• n+ii-) ■ B{<p -> e+e") = (3.14 ± 0.22 ± 0.14) • 10~ 4 was measured from the interfer- 
ence in the cross section of the process e + e~ — » The branching ratio B(<f> — > /i + /i~) = 
(3.30 ± 0.45 ± 0.32) • 10" 4 was obtained. 
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1 Introduction 



The <f> — > /i + fi~ decay reveals itself as an inter- 
ference pattern in the energy dependence of the 
cross section of process e + e~ — » in the 

region around ^-resonance peak. The interfer- 
ence amplitude is determined by the branch- 
ing ratio of the decay (j> — > The table 
value of the branching ratio B(<p — > /i + /i~) = 
(2.5 ±0.4) ■ 10~ 4 is based on the experiments 
on photoproduction of <f> meson ||, ||. In e + e~ 
collisions directly measurable is the square root 
of the product of the branching ratios 
yjB{^> -)• ' B((/) e+e~). The ^ - e- 
universality requires B(cj> — > /i + /j~) to be equal 
to B{4> — > e + e~) with the accuracy much higher 
than our experimental one despite the differ- 
ence in masses of the leptons. The measure- 
ments of the branching ratio B(4> — > e + e~) at 
e + e~ colliders are performed through summa- 
tion of cross sections over all decay channels 
|Q. The study of the interference e + e~ — > 
— ► (J. + gives the independent evaluation 
of the 0-meson leptonic width. First measure- 
ment of the decay <f> — > on e + e _ collider 
was performed in Orsay in 1972 ||. In this 
experiment the value of branching ratio of <j>- 
meson leptonic decay 

y/B(<f> fi+fi-) - B(<f> -»• e+e") = (2.93±0.96± 
0.32) • 10 -4 was obtained. Later similar mea- 
surements were performed in Novosibirsk || . 

2 Experiment 

The experiment was carried out with SND de- 
tector (Fig. at VEPP-2M in 1996-1997. SND 
is a general purpose non-magnetic detector [|| 
|{|. The main part of the SND is a spheri- 
cal electromagnetic calorimeter, consisting of 
1630 Nal(Tl) crystals. The solid angle of the 
calorimeter is ~ 90% of 4ir steradian. The an- 
gles of charged particles are measured by two 
cylindrical drift chambers covering 95% of Att 
steradian solid angle. The important part of 
the detector for the process under study is a 
muon system, consisting of streamer tubes and 
plastic scintillation counters, with 1 cm iron 
plates between blocks of tubes and counters. 
The simultaneous hits in the streamer tubes 
and scintillation counters produce a signal of 
muon system. The time difference between the 
hit and the beam collision is measured by inner 
JToj] and outer scintillation counters JTl| . 

The experiment was carried out in the en- 
ergy range 2E b = 984-1040 MeV and consisted 
of 6 data taking runs |lj, [ijj : 
PHL9601 - PHL9606. Five runs were used for 
analysis, corresponding to the total integrated 
luminosity AL — 2.61 pb _1 . 
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Figure 1: Detector SND — view across the beam; 
1 — beam pipe, 2 — drift chambers, 3 — in- 
ner scintillation counters, 4 — Nal(Tl) counters, 5 
— vacuum phototriodes, 6 — iron absorber, 7 — 
streamer tubes, 8 — outer scintillation counters 

About 4.6 • 10 6 of (f> mesons were produced. 
The integrated luminosity was measured with 
the accuracy of about 3% using e + e _ — » e + e~ 
and e + e~ — > 77 events. 

3 Event selection 

The energy behavior of the cross section of the 
process 

e+e- ^mV~(7) (1) 

was studied in the vicinity of <fi meson. Events 
with two collinear charged particles were se- 
lected for analysis. The selection of the events 
of the process (|l|) was performed with the fol- 
lowing cuts on angles of acollinearity of the 
charged particles in azimuth and polar direc- 
tions: I Acj> |< 10°, I A6 |< 25°. Additional 
photon emitted by either initial or final par- 
ticles was permitted. To avoid possible losses 
of events due to beam background or knock- 
on electrons in the drift chambers, one addi- 
tional charged particle was also permitted. To 
suppress the beam background the production 
point of charged particles was required to be 
within 0.5 cm from the interaction point in the 
azimuth plane and ±7.5 cm along the beam di- 
rection (the longitudinal size of the interaction 
region a z is about 2 cm). The polar angles of 
the charged particles were limited to the range 
45° < 9 < 135°, corresponding to the accep- 
tance angle of the muon system. 

The main sources of background are the 
cosmic muons and the processes 

e + e"-»e + e~, (2) 

e + e _ — > ir + ir~ (7). (3) 



To suppress the background from the process 
@ a procedure of e/ir(fi) separation was used. 
The algorithm is similar to that, developed for 
the ND detector It utilizes the difference 
in total energy depositions and the longitudi- 
nal energy deposition profiles for electrons, pi- 
ons and muons. As a result of this procedure 
the background from the process (||) was sup- 
pressed down to 4% of the events of the pro- 
cess (||). To suppress the contribution of the 
process ([|) a procedure of 7r//i separation by 
the muon system has been used. In the en- 
ergy range under study the probability to hit 
the muon system for a muon varies from 80 to 
93% and is as low as 1.5% for pions of the pro- 
cess ([}]). A requirement of a hit in the muon 
system also reduced the background from the 
process (^) by two order of magnitude, making 
its contribution negligible. 

After these cuts 80% of selected events are 
still cosmic background. The rejection of the 
events with the hit in two top segments of the 
muon system (about 45° in azimuth direction) 
(Fig. |l|) suppressed the cosmic events by two 
times. Further suppression of the cosmic back- 
ground was performed using the following pa- 
rameters: 




Figure 2: The time of the particle passage 
through the upper half of the muon system with 
respect to the beam collision time. 



1. the time difference between the hit of the 
muon system and the beam collision - t; 

2. the time difference between the hits in 
upper and lower halfs of the muon system 
- TOF; 

3. the sum of the distances from the tracks 
to the production point - R; 

4. the likelihood function , which is built 
on the basis of the energy depositions in 
the calorimeter layers for the muon: 

P fl = PA E l)- p A E 2)-P^(E 3 ), (4) 

where P lli (E i ) is the value of the proba- 
bility density function for the energy de- 
positions Ei in z-th calorimeter layer. These 
functions were obtained from the data 
sample, where muons were selected by 
muon system using strict cuts (| r |< 
5 ns, TOF < ns, see explanation be- 
low). 

The distribution of the parameter r is shown 
in Fig. 0. The scale for parameter r was ad- 
justed to have a peak at for e + e~ — > /i + 
events. The events with r |< 5 ns were se- 
lected for further analysis. In about 75% of 
the selected events there is a hit in the lower 
half of the muon system. For these events Fig. 
U shows the distribution of the time of flight 
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Figure 3: The time of flight TOF of the parti- 
cle from the upper half of the muon system to 
the lower one. The peak at is due to events of 
the e + e~ interaction, the right peak — cosmic 
background. 
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Figure 4: The RP distribution for the events 
with TOF < ns (hatched histogram) and 
events with TOF > 5 ns (empty histogram). 

TOF. Peaks in the spectra at TOF = ns 
and TOF = 7.2 ns are due to the process (0) 
and cosmic background respectively. 

To subtract the cosmic background the com- 
bination of the independent parameters R and 
Pfi was used: 

RP = ln(P p ) - 16 • R(cm) + 25. (5) 

The factor before R has been chosen to provide 
maximum separation between cosmic events and 
the events of the process ([!]). Additive con- 
stant equal to 25 fixes the scale of the param- 
eter RP. Fig. [| shows the RP distribution for 
the events with TOF > 5 ns (mainly cosmic 
events) and for the events with TOF < ns 
(mainly events of the process (0)). Most of 
the events of the process (|l|) are in the region 
RP > 0. At each energy point the number of 
events of the process (|l|) was estimated by the 
following formula: 



N n = N, 



RP>0 



C-N 



RP<0- 



(6) 



Here C = ffi F>0 ; N' RP>0 and iV^ p<0 were de- 
termined from the sample of the events with 
TOF > 5 ns; Nrp >0 , Nrp <0 — from the sam- 
ple of the events with TOF < 5 ns. 

The detection efficiency £ M for the process 
([!]) was determined from the Monte Carlo sim- 
ulation. The events generator was based on 
the formulae from the work |lq] . The passage 
of the particles through the detector was sim- 
ulated by the program UNIMOD2 |§. The 
energy dependence of the detection efficiency 
is especially important for the process under 
study. This dependence is determined mainly 
by the probability for a muon to reach the 
muon system. 



The detection efficiency determined from 
the Monte Carlo was corrected in order to take 
into account the effects, which were not in- 
cluded into simulation. First is the efficiency 
of a software 3-rd level trigger, which selected 
the events during the data taking. This effi- 
ciency equals to 97-98% for all runs, except 
the run PHL9601, where it was 93%. Second 
effect not included into the simulation is the 
inefficiency of the muon system. It was de- 
termined from the experimental data for each 
energy point and was also used as a correction 
for overall efficiency. This correction is about 
9% and was determined mainly by the dead 
time of the time channels. 



4 Data analysis 

The energy dependence of the detection cross 
section was fitted with the following formula: 

Ovis = E ll {E)a llll { ri ){E) + e w cr 7rTr ( 7 )( J B), (7) 

where E = 2E b ; e M (E) and e r « 0.017 — the 
detection efficiencies of the processes ((!]) and 
(H)i <t mm(7) an< ^ °V7r(7) — the cross sections of 
these processes. The energy dependence of the 
efficiency obtained from the simulation was ap- 
proximated by a smooth function. 

To perform a combined fit of the detec- 
tion cross sections for all experimental runs the 



scale factors e\ ^ for each run were introduced 
into the detection efficiencies as free fit param- 
eters. They are asymptotic detection efficien- 
cies at energies higher than 520 MeV, because 
the probability for a muon with such energy to 
hit the muon system is almost constant. This 
method provides an evaluation of efficiencies 
independent from simulation. 



The cross section a 



was calculated ac- 



cording to formulae from the work |L5| and tak- 
ing into account the experimental data on the 
pion form factor [|l7| and the decay <f> — > 2tt |jj . 
The cross section ct mai ( 7 ) (E) was taken in the 
following form: 



•W(7) 



(E) 



dE)- 13(E), 



(8) 



a^(E) = 83.50(n6)-f | Z 



Z = l-Q 



E 2 

TTIaTa 



E 2 - iET(E) ' 



where a flfl (E) is the Born cross section of the 
process e + e~ — > m^jT^ — mass and 

width of <p meson; Q, ip^ — module and phase 
of the interference amplitude; (3(E) — factor 
taking into account the radiative corrections. 
This factor was obtained as a ratio of the cross 



section of process (|l|) to the Born cross section 
for the same acceptance angles and with the 
0-meson contribution. The cross section of the 
process (Q) was calculated by the Monte Carlo 
method using the formulae from the work [[l5| 
with appropriate cuts on the angles and mo- 
menta of the final muons. 

The luminosity measurement was performed 
taking into account the interference term in the 
cross section of the process (||). The interfer- 
ence amplitude is about 0.5% and the phase 
differs by 180 degrees from the phase in the 
process ([!]). 

The fit of the experimental data gives the 
interference phase Vv = (4.5 ± 3.4)°, which is 
consistent with the theoretical expectation of 
0°. Therefore the final fit was made with the 
fixed phase ip^ = 0°. As a result the value of 
the interference amplitude 

Q = 0.129 ±0.009 

and the detection efficiencies e l s j for five runs 



were obtained (table [j]) . 



Table 1: 
process (II 



The detection efficiencies e l s j for the 



obtained from the simulation and 
the experiment. 



No 



£\f,MC 



£ 



f , exp 



PHL9601 
PHL9602 
PHL9604 
PHL9605 
PHL9606 



0.270±0.003 
0.282±0.004 
0.287±0.003 
0.289±0.003 
0.287±0.003 



0.254±0.005 
0.277±0.004 
0.292±0.003 
0.290±0.003 
0.260±0.003 



One can see from the table that the detec- 
tion efficiencies, obtained from the simulation 
and the experiment, are in good agreement, ex- 
cept for the runs PHL9601 and PHL9606. In 
PHL9601 run the muon system worked in dif- 
ferent regime. It was not fully compensated by 
the corrections. In addition, the efficiency of 
the software trigger in this run was lower. The 
explanation of the lower efficiency in the run 
PHL9606 is a lower gain in the drift chambers 
and consequently a worse resolution in 9. 

The error in the detection efficiency does 
not directly contribute into the error at inter- 
ference amplitude if each run is recorded in the 
same conditions because the interference am- 
plitude is a relative value. Therefore the data 
with low gain in the drift chambers were ex- 
cluded from the runs. 

Fig. H shows the Born cross section of the 
process e + e~ — * [i + [i~ . The x 2 value obtained 
in the fit equals to 36.0 for 45 degrees of free- 
dom. 



The systematic error of the amplitude of 
the interference is determined by the errors of 
the luminosity estimation and the calculation 
of the radiative corrections. The measurement 
of the luminosity using the events of the pro- 
cess e + e~ — > 77 gives the estimation of the sys- 



tematic error Sh 



2.1%. Additional system- 



atic error can be due to imprecise evaluation of 
the background from the process (||) . The esti- 
mation of this contribution is = 0.6%. The 
contribution of the error in the calculation of 
the radiative corrections to the systematic er- 
rors of the interference amplitude is Sp = 4%. 
The resulting systematic error is S = 4.6%. 
The interference amplitude is related to the 



branching ratio of the decay 
following formula: 



by the 



Q = 



3^B( 



-)B{ 



(9) 



From 



where a is the fine structure constant, 
this relation the product is 

(3.14±0.22± 



0.14) 

e+ 
B( 



1-4 



10~ 4 . Using the table value £?(</> — > 
(2.99 ± 0.08) • 10~ 4 [@, we obtain 
= (3.30 ±0.45 ±0.32) • 10~ 4 . 



5 Conclusion 

In this work the energy dependence of the cross 
section of the process e + e~ — > in the 

vicinity of the <p resonance was studied. The 
interference pattern determined by the decay 
— > /J + /i~ was measured, giving the product 
of the leptonic branching ratios of <j> meson 
y/B{<j> -> e+e-)B((j> = (3.14±0.22± 

0.14) • 10~ 4 . Assuming ^ — e-universality one 
can compare this value with the table branch- 
ing ratio B((f> -> e+e~) = (2.99 ± 0.08) • 10~ 4 
|l[ They are in good agreement. Not using 
/i — e-universality the branching ratio B(cj> — > 
H + = (3.30 ± 0.45 ± 0.32) • 10~ 4 was ob- 
tained. It is close to one standard deviation 
from the table value B(<f> — > = (2.5 ± 

0.4)-10- 4 §. 
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